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ABSTRACT 


The  historical  development  of  the  foundations  of  thermomechanics  is 
reviewed,  with  emphasis  placed  upon  the  critical  contribution  of  each  in¬ 
vestigator.  The  theories  are  discussed  in  terms  of  their  relevance  to 
real  materials  in  the  important  high-velocity  regime  which'  lies  between 
the  low  and  intermediate  velocity  regime  (which  involves  essentially  iso¬ 
thermal  elastic  and  plastic  wave  propagation)  on  one  hand  and  the  ultrahig’n- 
velocity  regime  (in  which  a  hydrodynamic  description  of  the  material  be¬ 
havior  is  adequate)  on  the  other.  Various  facets  of  high-velocity  wave 
propagation  in  various  materials  involve  elastic,  plastic,  viscous,  and 
thermal  effects  on  a  macroscopic  scale  and  perhaps  on  a  microscopic  scale. 
Evaluations  are  made  from  the  standpoint  of  thermodynamics  as  well  as 
classical  mechanics. 

Using  fundamental  principles  of  thermodynamics  and  mechanics,  thermo- 
mechanically  coupled  linear  energy  and  constitutive  equations  are  derived 
for  the  elastic  solid  with  microunit  cells  discussed  by  Mindlin. 

Equations  of  motion  for  the  microstructure  are  found  for  the  large 
amplitude  case  where  Eulerian  and  Lagrangian  formulations  are  not 
identical. 
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I.  REVIEW  OF  PREVIOUS  WORK 

A.  INTRODUCTION 

In  the  past  decade,  a  very  large  body  of  experimental  data  on  the 

behavior  of  various  solids  under  strong  shock  loading  has  been  obtained 

12  3 

by  means  of  a  variety  of  high-pressure  dynamic  loading  techniques.  *  ’ 

Since  most  of  the  original  experimental  work  was  conducted  at  very  high 
pressures  (say,  three  or  more  orders  of  magnitude  above  the  ultimate 
strength)  at  which  it  has  been  recognized  that  strength  effects  are  not 
important,  the  analyses  developed  to  describe  this  work  mathematically 
were  based  on  hydrodynamics  of  dense,  inviscid  fluids.  More  recent  effort 
has  been  devoted  to  hypervelocity  impact  (HVI) ,  i.e.,  the  impact  of  a  small 
body  with  a  massive  body  at  relative  velocities  exceeding  the  acoustic 
velocities  cf  both.  Since  HVI  involves  creation  of  craters  of  finite  size, 
it  must  depend  upon  the  strength  properties  of  the  materials  involved. 

Most  of  the  attempts  to  analyze  HVI  have  involved  either  minor  modification 
of  hydrodynamic  theory  or  correlation  of  highly  simplified  phenomenological 
models  with  experimental  results. 

An  example  of  the  inadequacy  of  the  hydrodynamic  approach  for  the 

4 

high-velocity  regime  has  been  reported  by  Curran.  The  results  of  his 
plate-slap  experiments  on  2024  aluminum-alloy  plates  at  an  impact  velocity 
of  1.9  ram/ps  indicated  that  attenuation  begins  approximately  80-percent 
sooner  and  is  approximately  twice  as  severe  as  predicted  by  hydrodynamic 
theory.  An  analysis  able  to  predict  this  greater  attenuation  accurately 
is  highly  desired. 

Until  very  recently,  thermodynamic  aspects  of  high-velocity  wa/e 
propagation  have  been  highly  ne.’lected.  Although  the  theory  of  linear 
elastic  materials  with  thermomechanical  coupling  (through  the  energy 
equation)  dates  back  to  1837  with  Duharamel's  work,^  there  have  been  few 
new  major  accomplishments  in  thermomechanics  until  about  ten  years  ago. 
Since  the  topic  of  linear  thennoelasticity  without  thermomechanical 

6 

coupling  has  been  treated  extensively  in  Boley  and  Weiner's  textbook, 
it  is  omitted  here. 
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B.  THEPMOELASTICITY 

Although  the  energy  equation  and  constitutive  equations  for  a 
classical  linear  elastic  solid  with  thermomechanical  coupling  were 
originated  by  Duhammel  in  1837,  surprisingly  little  work  has  been  done 
in  this  area.  Here  the  term  "classical"  means  one  without  micros tructural 
effects  and  with  elastic  properties  independent  of  temperature.  Most  of 
the  analyses  reported  in  the  literature  are  concerned  with  the  response 
(thermal  and  mechanical)  due  to  sudden  or  cyclic  heating  but  no  mechanical 
loading.  It  appears  that  more  analysis  should  be  devoted  to  problems 
involving  purely  mechanical  loading,  with  particular  emphasis  on  determin¬ 
ing  the  loading  rate  at  which  the  thermomechanical  coupling  effects 
become  important.  On  a  strictly  intuitive  basis,  it  is  believed  that  in 
the  regime  in  which  thermomechanical  coupling  effects  become  important, 
the  use  of  a  finite  description  of  strain,  rather  than  the  linear 
infinitesimal  one,  is  necessary. 

One  of  the  most  interesting  anomalies  in  the  behavior  of  a  classical, 
linear  elastic  solid  with  thermomechanical  coupling  has  been  discussed  by 
Deresiewicz:^  His  analysis  predicted  that  the  phase  velocity  of  plane 
elastic  waves  approach  the  adiabatic  value  at  low  frequencies  and  the  iso¬ 
thermal  value  at  high  frequencies.  Deresiewicz  pointed  out  that  this  is 
in  disagreement  with  the  Laplace-Rayleigh  thermomechanical  theory  of 
gases,  which  Deresiewicz  claimed  is  erroneous. 

8—10 

Although  various  theories  of  elasticity,  taking  into  considera¬ 
tion  doublet  stresses  (doublet  forces/area)  in  addition  to  classical 
stresses  (forces/area),  have  been  proposed  to  explain  such  phenomena  as 
size  effects  in  materials  and  crystal-lattice-type  wave  dispersion;  only 
recently  has  attention  been  given  to  the  thermodynamic  aspects  of  such 

the  theories  of  thermoelasticity  known  to  the  authors, 
dependence  of  the  elastic,  cons  ants  is  neglected,  even 
thermomechanical  theory  of  elastoplastic  waves  in  solids 


solids 
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In  all  of 
the  temperature 
though  a  recent 
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considers  an  effect  of  the  same  order  of  magnitude:  the  effect  of  tempera¬ 
ture  on  the  yield  strength. ^  Rosenfield  and  Averbach, ^  in  an  analysis 
which  is  on  a  rather  weak  mathematical  and  thermodynamic  basis,  indicated 
that  the  temperature  dependence  of  the  elastic  coefficients  is  related  to 
the  stress  dependence  of  the  thermal-expansion  coefficient  a  as  follows: 

a  =  a  +  a,  a  (1) 

o  1 

where 

a-  -  -E-2  dE/dT  (2) 

where  aQ  is  the  coefficient  of  thermal  (linear)  expansion  at  zero  stress, 
is  the  change  in  this  coefficient  per  unit  of  normal  stress,  o  is 
the  normal  stress  applied  in  the  direction  in  which  a  is  measured,  E 
is  Young's  modulus,  and  T  is  temperature.  Rosenfield  and  Averbach  con¬ 
ducted  experiments  which  confirmed  the  linear  dependence  of  a  with 
stress,  as  predicted  by  Eq.  1.  However,  the  discrepancy  between  tne  value 
of  the  proportionality  constant  predicted  by  their  theory,  Eq.  2,  and 

that  measured  varied  from  6  to  100  percent.  An  analysis  which  is  on  a 
more  sound  thermodynamic  and  mathematical  basis  than  that  of  Ref.  16  is 
given  in  the  appendix.  Unfortunately  the  new  theory  requires  additional 
experimental  data  which  are  not  normally  available,  namely  the  variation 
of  with  temperature. 

6  17 

Following  a  procedure  suggested  in  Boley  and  Weiner's,  Fillon 
formulated  a  theory  of  thermoelasticity  for  physically  nonlinear  elastic 
solids.  He  considered  the  strains  as  independent  thermodynamic  variables 
("fluxes")  and  the  stresses  as  dependent  thermodynamic  variables  ("forces"); 
thus,  he  utilized  the  Helmholtz  free  energy  function  (also  called  the 
Gibbs  work  function).  It  appears  that  he  might  just  as  well  have 
considered  the  stresses  as  independent  variables  and  the  strains  as 

18 

dependent  variables,  and  then  utilized  the  Gibbs  thermodynamic  potential. 
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Vakulenko  has  considered  stress  as  the  independent  thermodynamic  variable 

20 

in  his  work  in  thermoplasticity,  and  Burridge  and  Knopoff  have  introduced 
a  formulation  of  mechanical  energy,  using  stress  as  the  independent 
variable,  which  is  valid  for  prestressed  and  inelastic  media  as  well  as 
nonp res  tressed  elastic  media.  This  approach  is  used  in  the  appendix  to 
considerably  simplify  the  derivation  of  the  stress  dependency  of  the 
thermal-expansion  coefficient. 

Another  type  of  nonlinearity,  namely  geometric  nonlinearity  (i.e., 

21 

finite  deformations),  has  been  discussed  very  ably  by  Lee  and  Liu,  and 

some  of  the  underlying  thermodynamics  aspects  have  been  put  on  a  sound 

22  23 

mathematical  basis  by  Coleman  and  Noll  and  Wang  and  Bowen. 

C.  THERMOPLASTICITY 

24 

The  quantitative  experiemnts  of  Farren  and  Taylor  in  1925  form  the 
basis  for  thermoplasticity.  They  found  that  approximately  90  percent  of 
the  work  done  in  plastic  deformation  is  dissipated  in  the  form  of  heat  and 
the  remaining  10  percent  is  irreversibly  stored  due  to  permanent  distortion 
of  the  crystal  lattice. 

Some  of  the  work  in  thermoplasticity  has  been  largely  of  a  phlloso- 
25-27 

phical  nature.  The  first  attempts  to  put  the  subject  on  a  quantitative 

28  29  30 

basis  were  made  independently  by  Freudenthal,  Prager,  Ziegler,  and 

31  19  32 

Vakulenko,  *  Their  work  has  been  discussed  critically  by  Naghdi, 

who  considered  their  work  strictly  exploratory  in  nature.  His  primary 
objection  was  that  their  theories  are  uncoupled,  i.e.,  they  do  not  con¬ 
tain  a  heat  conduction  term.  He  also  questioned  their  choice  of  state 
variables  and  generalization  of  Onsager's  principle  to  nonlinear 
phenomena. 

33 

Grigorian  took  the  viewpoint  that,  although  thermodynamic  concepts 
are  of  limited  usefulness  in  constructing  the  mechanical  parts  of  consti¬ 
tutive  equations,  one  ohould  check  the  thermodynamic  consequences  of  a 
tentative  constitutive  relation  in  order  to  make  the  mathematical  model 
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thermodynamically  complete.  (This  point  is  illustrated  in  the  consideration 
in  the  appendix.)  It  is  interesting  to  note  that  Grigorian  recognized 
that  some  of  the  plastic  work  goes  into  increasing  the  internal  energy  due 
to  changes  in  structure  of  the  material. 

34 

Foster  devoted  attention  to  relaxation-spectrum  effects.  In 
particular,  he  noted  that  there  are  different  relaxation  spectra  for 
loading  (increasing  load)  and  unloading  (decreasing  load) . 

Assuming  that  all  of  the  plastic  work  is  dissipated  in  the  form  of 

35 

heat,  Dillon  formulated  a  theory  of  thermoplasticity  with  thermomechanical 

36 

coupling  for  small  strains  and  Backman  formulated  an  analogous  theory 

for  finite  strains.  However,  the  latter  work  has  been  criticized  and 

21 

improved  upon  by  Lee  and  Liu. 

To  date  the  most  complete  mathematical  models  for  dynamic  thermo- 

37 

plasticity  analyses  are  those  of  Grigorian  (with  emphasis  on  soils)  and 
Lee  and  Wierzbicki^  (with  particular  reference  to  metals). 

As  mentioned  previously,  Lee  and  Wierzbicki^  included  the  effect 

of  temperature  on  the  yield  strength,  but  neglected  the  thermodynamically 

related  effect  of  stress  on  the  thermal-expansion  coefficient,  a  .  In 

16 

their  experiments,  Rosenfield  and  Averbach  measured  this  effect  for 
stresses  well  into  the  plastic  range.  They  found  that  as  soon  as  the 
elastic  limit  is  exceeded,  a  increases  with  stress  at  a  greater  slope 
an  in  the  elastic  range.  However,  the  slope  decreases  with  increasing 
stress,  eventually  becoming  zero,  and  finally  negative.  When  they 
removed  the  stress  after  each  successive  loading,  they  found  that  the 
permanent  changes  in  a  also  increased,  then  leveled  off,  and  finally 
decreased  at  higher  applied-stress  levels.  This  appears  to  be  closely 
related  to  the  observations  of  Farrsn  and  Taylor  that  the  percentages 
oi  the  total  plastic  work  which  goes  into  permanent  distortion  of  the 
mice  varies  with  applied  strain.  In  view  of  these  considerations,  it 
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is  conjectured  that  the  effect  of  the  stress  dependency  of  the  thermal- 
expansion  coefficient  would  be  greatest  for  short-duration  pulses  in 
which  the  loading  is  removed  before  there  is  sufficient  time  for  very 
large  plastic  strains  to  be  reached.  However,  as  in  the  case  of  the  same 
effect  in  the  elastic  range,  quantitative  inclusion  of  this  effect  in  a 
theory  must  await  experimental  measurements  of  it  over  a  range  of 
temperatures  ranging  upward  from  room  temperature. 


D.  THERMOVISCOELASTICITY 

The  theory  of  thermoviscoelasticity  has  strong  thermodynamic  bases, 

38-40  41-42  43-45 

thanks  to  extensive  research  by  Biot,  Chu,  Schapery,  and 

Coleman. ^6-/t8  other  research  on  the  formulation  of  constitutive  equations 

49  50 

for  thermoviscoelasticity  has  been  carried  out  by  Eringen,  Cowin, 

51  52 

Chudnovskii,  and  Leigh. 


Various  aspects  of  the  nonlinear  theory  of  viscoelasticity  have 

been  treated  in  recent  work  by  Koh  and  Eringen, Kline, ^  and  Valanis.^ 
cc 

Recently  Eringen^'  has  formulated  a  theory  of  micros tructural  viscoelasticity 
including  thermodynamic  considerations. 

The  literature  which  has  appeared  in  the  past  decade  on  dynamic 
thermoviscoelastic  (DTVE)  problems  is  too  voluminous  to  mention  in  detail 
here.  For  a  basic  treatment  of  the  subject,  reference  is  made  to  Boley 
and  Weiner's  book.6  Most  of  the  DTVE  analyses  found  in  the  literature 
neglect  thermomechanical  coupling.  An  example  of  one  of  the  few  DTVE 
analyses  in  which  this  coupling  is  considered  is  found  in  Ref.  57. 

It  appears  that  the  possibility  cf  the  type  of  thermomechanical 

CO 

instability  discussed  by  Berg  should  be  investigated  further,  especially 
for  high-intensity  waves.  He  pointed  out  that  the  softening  of  a  body 
due  to  heating  allows  external  loads  to  deliver  more  power  to  the  body  at 
an  increasing  rate. 
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E.  WAVE  PROPAGATION 

There  are  two  fundamentally  different  types  of  waves: 

1.  Simple,  continuous  waves,  which  can  be  mathematically  described 
by  a  Fourier  series  of  harmonic  waves.  These  waves  are  of 
small  amplitude  and  are  often  called  acoustic  waves. 

2.  Discontinuous  waves  which  involve  a  discontinuous  jump  in 

some  physical  quantity.  The  order,  n  ,  of  a  discontinuous 

wave  is  the  order  of  the  time  derivative  of  the  displacement 

which  is  discontinuous.  Thus,  a  shock  wave,  which  involves 

a  jump  in  velocity,  is  a  first-order  wave  and  is  usually  the 
* 

most  severe.  Higher  order  waves  involve  jumps  in  accelera¬ 
tion  (acceleration  waves)  and  in  higher  time  derivatives. 


59 

The  classical  treatment,  as  used  by  Rayleigh,  is  adequate  for 
handling  acoustic  waves.  The  literature  on  the  propagation  of  such  waves 
in  elastic  and  viscoelastic  solids  is  too  voluminous  to  mention  here 

except  to  refer  to  two  texts60,61  and  three  review  papers  devoted  to  this 

.  4  £2-64 

topic. 


The  common  metallic  alloys  used  in  load-carrying  applications, 
namely  steel,  aluminum  and  most  titanium  alloys,  behave  macros copically 
as  isotropic  materials  at  least  in  the  elastic  range.  However,  since 
they  are  really  aggregates  of  many  highly  anisotropic  crystals,  they  behave 
anisotropically  on  a  microscale.  This  has  been  the  motivation  for  various 

65_68  ^9 

theoretical  and  experimental  studies  of  acoustic-wave  propagation 

in  single  crystals.  A  few  metals  and  alloys  of  increasing  engineering 
importance,  namely  beryllium,  certain  titanium  alloys,  zirconium,  and 
uranium,  are  anisotropi.:  even  on  a  macroscale.  Furthermore,  filamentary 
composites  behave  anisotropically  on  a  macroscale. 

ft 

An  alternative,  but  equivalent,  definition  of  the  order  n  of  a  discon¬ 
tinuity  is  one  plus  the  order  of  the  space  derivative  of  the  strain  which 
is  discontinuous.  Thus,  a  shock  wave,  which  involves  a  jump  in  strain, 
i,  a  first-order  wave. 
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The  development  of  theories  of  solids  with  microstructure^*^’^ 
now  permits  a  mathematical  treatment  of  solids  which  are  isotropic  on  a 
macroscale  and  yet  anisotropic  on  a  microscale.  This  opens  a  new  avenue 
of  improved  correlation  between  theory  and  experimental  studies  of 
acoustic  wave  propagation.  However,  to  date,  studies  along  this  line 

10  70-72 

have  been  limited  to  theoretical  studies  of  relatively  simple  systems.  ’ 

The  method  of  characteristics,  which  is  based  on  the  concept  of  a 

characteristic  surface  or  discontinuity  surface  appears  to  have  been 

73 

developed  at  the  turn  of  the  century  by  Hadamard  for  the  study  of  the 
propagation  of  discontinuities  in  hydrodynamics.  The  method  has  found 
wide  use  in  fluid  mechanics ,  including  gas  dynamics  and  magnetofluidmech- 
anics  as  well  as  hydrodynamics . 

Apparently  the  first  application  of  the  method  of  characteristics 

to  solids  behaving  nonhydrodynamically  is  due  to  Thomas.  He  was  originally 

74 

concerned  with  plastic  solids  but  later  generalized  his  compatibility 
75  76  77 

conditions.  Hill//  discussed  stability  and  uniqueness  of  acceleration 

waves  in  general  media  with  emphasis  on  incompressible  and  elastoplastic 

media.  Acceleration  waves  in  anisotropic  elastic  solids  were  analyzed  by 

78  79 

Howard,  who  was  concerned  with  transversely  isotropic  media  and  Nariboli 

who  treated  general  anisotropic  elastic  media.  Nariboli's  approach  is 
interesting  in  that  he  combined  the  singular  surface  concept  of  Thomas 
with  the  ray  theory  of  Courant  and  Hilbert  to  obtain  an  equation  which  up¬ 
on  integration  predicts  the  growth  or  decay  of  the  discontinuities.  Very 
80 

recently  Lur'e  has  analyzed  the  propagation  of  discontinuities  in 
continuation  with  microstructure. 

Recently  a  number  of  investigators^  ^  have  investigated  nonlinear 
aspects  of  acceleration  waves  and  higher-order  waves  in  nonlinear  elastic 
materials.  Some  of  the  general  results  of  these  studies  are: 

1.  When  the  strain  has  a  discontinuity,  the  derivative  of  the 

82 

strain  also  must  have  a  discontinuity.  Thus,  shock  waves 
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(to  be  discussed  later)  are  accompanied  by  acceleration  waves 
although  the  converse  situation  does  not  necessarily  hold. 

83 

2.  Varley  and  Dunwoody  found  that  the  linear  theory  of  irrational 

waves  has  an  error  which  increases  with  time.  Their  nonlinear 

theory  predicts  that  either  a  shock  forms  or  the  material 

"forgets"  the  details  of  the  disturbance  in  a  finite  time. 

84 

Green's  theory  predicts  that  acceleration  waves  either 
become  infinite  in  finite  time  (i.e.,  a  shock  forms)  or  decay 
to  zero  in  infinite  time,  with  the  exception  of  transverse 
waves  normal  to  a  shock. 

3.  Discontinuities  of  higher  order  than  the  second  propagate 

84 

with  constant  strength  and  at  the  same  velocity  as  the 

,  .  85 

acceleration  wave. 

86 

In  their  book,  Truesdell  and  Noll  presented  a  review  of  accelera¬ 
tion  waves  in  generalized  elastic  media.  Recently  considerable  attention 

87-93 

has  been  devoted  to  acceleration  waves  in  viscoelastic  media.  Of 

these,  the  work  of  Varley  is  especially  important  because  the  material  may 
be  generally  anisotropic  and  nonhomogeneous . 

The  problem  of  the  propagation  of  discontinuities  into  a  prestressed 

medium  is  important  in  connection  with  propagation  of  unloading  waves  (if 

the  duration  of  loading  is  sufficiently  long  compared  to  the  transit  time 

of  the  reflected  wave),  since  then  the  medium  is  essentially  prestressed 

by  the  loading  wave.  For  small  deformations  (i.e.,  acoustic  waves),  this 

94 

problem  was  treated  by  Biot.  The  case  of  acoustic  waves  but  large 

deformation  due  to  the  prestress  have  been  considered  in  Ref.  95-97. 

Finally,  the  propagation  of  finite  discontinuities  in  prestressed  media 

98 

has  been  treated  by  Green. 

Although  the  foregoing  analyses  employ  thermodynamics  in  the 
formulation  of  the  constitutive  relations,  they  do  not  consider  thermo¬ 
mechanical  coupling.  This  has  been  treated  in  two  recent  papers  by 
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Chadwick  and  Powdrill  for  the  linear  thermoelastic  case  and  by  Fine 

for  the  thermoelastoplastic  case. 

In  the  hydrodynamic  formulation  of  the  shock  wave  problem,  the 
major  physical  effect  considered  in  the  constitutive  equation  is 
compressibility,  i.e.,  dilatational  (volumetric)  waves.  Thus,  shear  waves, 
in  fact  any  shear  effects  whatsoever,  are  not  considered  at  all.  As  has 
been  pointed  out  by  Fyfe  et  al,^^  this  does  not  mean  that  the  hydrodynamic 
theory  cannot  predict  some  physical  aspects  of  a  shock  wave  even  at  high 
velocities  (as  opposed  to  hypervelocities).  An  example  of  an  early  attempt 
to  combine  shear  effects  with  dilatational  effects  is  contained  in  Ref.  102. 

It  should  be  pointed  out  that  shock  waves  are  not  limited  to 

compressible  media.  In  fact,  recent  analyses  by  Chu^^  ^  and  by 
105 

Collins  have  been  concerned  with  the  propagation  of  shock  waves  in 

* 

incompressible  isotropic  elastic  media.  In  such  media,  only  nondilatational 
shock  waves  can  occur. 

Recently  Bland, Green, ^  and  Coleman,  Gurtin,  and  Herrera, ^ 

considered  plane  shock  waves  in  compressible  isotropic  elastic  media. 

Although  these  investigators  considered  certain  thermodynamic  aspects, 

108 

they  did  not  consider  them  to  the  depth  considered  by  Dewey  and  by 

Lee  and  Wierzbicki.^  The  case  of  plane  shock  waves  in  compressible 

isotropic  thermoelastic  media  was  recently  treated  by  Chadwick  and 
109 

Powdrill.  Their  analysis  is  interesting  in  that  they  claim  that  their 
analysis  leads  to  an  infinite  family  of  shock  wave  solutions  which  range 
thermodynamically  from  isothermal  to  adiabatic.  They  claim  that  they  know 
of  no  thermodynamic  or  mechanical  principle  which  enables  selection  of 
the  "correct"  solution,  although  they  suggest  the  possibility  of  using 
the  theory  of  nonlinear  hyperbolic  partial  differential  equations  or 
the  irreversibility  principle  of  thermodynamics  (second  law).  The 
present  writers  suggest  that  a  stability  investigation  might  be  in  order 
here,  i.e.,  it  is  unlikely  that  all  of  the  possible  solutions  are  stable. 

A 

For  an  isotropic  linearly  elastic  solid,  incompressibility  requires 
that  Poisson's  ratio  be  equal  to  1/2.  This  value  is  closely  approached 
in  certain  rubber-type  materials,  but  not  in  metallic  materials. 
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In  Ref.  106,  Bland  found  that  two  kinds  of  shock  waves,  dilatational 
and  nondilatational,  are  possible  in  a  compressible  isotropic  elastic 
solid.  However,  it  is  interesting  to  note  that  he  concluded  that  nondilata¬ 
tional  shock  waves  cannot  propagate  into  a  compressible  hyperelastic  solid 
at  rest  in  its  reference  state.  To  the  present  writers,  this  implies  that 
nondilatational  shock  waves  can  occur  only  during  unloading,  since  then 
the  material,  being  already  loaded,  is  not  at  its  reference  state.  In 
Ref.  107,  Bland  showed  that  in  many  solids  (including  simple  elastic  mate¬ 
rials),  the  only  kind  of  shock  waves  which  are  stable  and  admissible 
thermodynamically  are  tensile  dilatational  shock  waves  (in  contrast  to 
compressive  dilatational  shock  waves  in  fluids).  Although  plane  compres¬ 
sive  shock  waves  in  elastoplastic  media  have  been  studied  extensively, 

Bland  claimed  that  they  are  possible  only  in  media  in  which  the  stress- 
strain  curve  is  concave  upward  in  the  plastic  range. 

To  study  the  shock  structure  (i.e.  the  shape  of  the  shock  wave 
front),  it  is  necessary  to  include  viscous  dissipative  effects.  Recently 
Bland‘S  has  performed  such  an  analysis  for  plane  waves  in  a  simple 
viscoelastic  solid.  He  showed  that  all  monotonic  wave  profiles  ultimately 
adopt  the  same  constant  wave  profile,  and  he  calculated  the  profile  width 
and  the  time  formation  of  the  profile. 


In  this  brief  review,  a  number  of  related  special  topics  should  at 
least  be  mentioned: 

1.  Disintegration  or  pulverization  of  material  behind  a  shock 
front,  applicable  to  certain  soils. 

112 

2.  Changes  of  phase  such  as  melting  due  to  shock  wave  heating, 

113 

shock-induced  metallurgical  transformations,  and  shock 
hardening  of  metals. 


3.  Interactions  between  acoustic  waves  and  shock  waves, 


115 


reflection  and  refraction  of  shock  waves  at  free  surfaces  and 
interfaces,  and  spallation  fracture  criteria. 


Experimental  techniques  for  generating  shock  waves  and  measur- 

1-3 

ing  various  characteristics  of  them. 
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II. 


COUPLED  THEKMOELASTICITY  THEORY  FOR  A  LINEAR  SOLID  WITH 


MICROSTRUCTURE 

A.  INTRODUCTION 

The  use  of  high-pressure  dynamic  loading  techniques  in  studying  the 

behavior  of  solids  under  strong  shock  loading  has  produced  a  very  large 

3 

body  of  experimental  data  in  recent  years.  Most  of  the  original 
experimental  work  was  associated  with  very  high  pressures  (three  or  more 
orders  of  magnitude  above  the  ultimate  strength  of  the  material)  at  which 
it  has  been  recognized  that  strength  effects  are  not  important.  Thus, 
the  theory  developed  to  accompany  this  work  was  based  on  hydrodynamics 
of  dense  fluids.  Later  interest  turned  to  hypervelocity  impact  of  a  swell 
body  on  a  massive  body.  This  problem  involved  creation  of  craters  of 
finite  size  and  consequently  involved  high-pressure  waves  that  must 
necessarily  depend  on  the  strength  properties  of  the  materials  involved. 
The  natural  extension  of  hydrodynamic  theory  to  this  problem  has  been 
by  the  most  expedient  addition  of  special  effects  to  existing  computer 
codes  (programs)  and  in  correlation  of  simple  phenomenological  models 
with  experimental  results. 

There  is  a  well  recognized  need  for  an  attack  on  the  problem  of 
propagation  of  high-pressure  waves  form  a  fundamental  thermomechanical 
point  of  view.15  This  effort  is  underway  in  several  areas.  The  effort 
reported  here  is  concerned  with  developing  an  applicable  coupled 
thermoelastic  theory  that  will  take  into  consideration  the  microstructure 
of  an  elastic  solid. 


The  use  of  "theory"  here  is  in  its  most  general  sense,  since  most  of 
what  has  been  done  in  the  past  is  in  the  nature  of  computer  solutions 
to  very  involved  nonlinear  flow  problems  and  might  properly  be  called 
computer  experiments  on  mathematical  simulations. 
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The  coupled  thermomechanical  theory  (thermal  and  mechanical  effects 
coupled  in  energy  and  constitutive  equations)  is  well  established  for  the 

macroscopic  description  of  solid  media.  The  elastic  case  is  due  to 

5  AO 

Duhammel  in  1837,  the  viscoelastic  case  to  Biot  in  1955,  and  the 

35 

viscoplastic  case  to  Dillon  in  1963. 

Experimental  observation  of  such  mechanical  phenomena  as  size 
effects  and  crystal-lattice  type  wave  dispersion  effects  has  led  to  various 
micromechanical  theories.  Of  these,  perhaps  the  most  notable  is  the 
couple-stress  theory,  which  in  addition  to  ordinary  stresses  (forces  per 
unit  area)  also  considers  couples  per  unit  area.  This  theory  was 

g 

originated  by  Voigt  in  1887  and  was  considerably  developed  by  the 

9 

Cosserat  brothers  approximately  sixty  years  ago.  In  the  past  ten  years, 
this  theory  has  been  developed  further  by  various  investigators  too 
numerous  to  mention  here. 

Mindlin^  noted  that  couple  stresses  are  actually  only  the  antisymet- 
ric  part  of  a  more  general  tensor  which  he  called  the  double  stress  tensor; 
thus,  he  originated  a  theory  which  he  called  linear  elasticity  with  micro¬ 
structure.  He  also  pointed  out  that  hie  theory  is  mathematically 

equivalent  to  a  linearized  version  of  Ericksen  and  Truesdell's  theory  of 

118 

deformable  directors.  In  the  latter  theory,  the  mechanical  behavior 
of  a  body  at  a  given  point  is  assumed  to  depend  upon  not  only  the  deforma¬ 
tion  of  the  point  but  also  the  deformaticn  of  an  oriented  "director" 
located  at  the  point.  If  the  director  is  assumed  to  be  rigid,  the 
equations  reduce  to  those  of  the  Voigt-Cosserat  couple-stress  theory. 

Upon  completion  of  the  present  work,  the  work  of  Green,  Rivlin,  and 
119-121 

Naghdi  came  to  the  attention  of  the  authors.  In  Refs.  119  and  121 

they  derived  an  energy  equation  for  a  simple  multipolar  material,  but 
they  did  not  deal  with  constitutive  equations.  In  Ref.  120  they  extended 
their  previous  work  to  include  multipolar  deformation  fields.  Their  work 
differs  from  ours  in  that  they  started  with  an  entirely  different  set  of 
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hypotheses  in  their  theory.  We  believe  that  our  derivation  on  the  basis 
of  a  microelement  is  more  satisfactory  from  a  physical  viewpoint.  Also, 
they  did  not  include  initial  stresses  nor  anisotropic  behavior  in  their 
constitutive  equation. 

In  the  work  reported  here,  the  usual  material  constitutive  equations 
are  replaced  by  equations  derived  for  the  specific  case  of  shock  wave 
loading.  This  is  consistent  with  the  current  interest  in  constitutive 
equations  that  have  been  experimentally  determined  from  plate  slap  and 
other  strong  shock  wave  techniques. 

B.  CONSERVATION  OF  ENERGY 

Using  Mindlin's  expression"^  for  the  kinetic  energy  density,  the 
kinetic  energy  T  of  a  macrovolume  V  can  be  rewritten  as  follows  , 


T 


A  /  pijij  dv  +  A  JV  dv 

V  V 


(3a) 


where  p  =  +  p ' 


(3b) 


d  U 
q\ 


pl6ql4Al*jl 


6p26q2^2&j2 


p3  q3  £3  j3 


■) 


(3c) 


where 


displacement  components 

semilengths  of  edges  of  microelement 

microdeformation 

mass  of  macromaterial  per  unit  macrovolume 
mass  of  micromaterial  per  unit  macrovolume 
K rone eke r  delta 

direction  cosines  of  edges  of  microelement  with  respect 
to  fixed  spatial  coordinates 
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and  a  dot  denotes  differentiation  with  respect  to  time.  The  first  term 
on  the  right-hand  side  of  Eq.  (3a)  is  the  classical  kinetic  energy 
expression  and  the  second  term  represents  the  contribution  of  the  micro¬ 
element. 

Mindlin's^  expression  for  the  rate  of  work  W  done  on  the  macro¬ 
volume  V  by  the  external  forces  (surface  and  body  types)  is  equivalent 
to  the  form: 


W  -  J dA  +  J fjij  dV  ♦  Jy 


* jk  dA  + 


/VjkdV  « 


where  t^  =  surface  force  per  unit  area 

f  =  body  force  per  unit  macrovolume 

T..  =  surface  force  doublet  per  unit  area  with  moment  arm 
jk 

normal  to  axis  X.  and  acting  in  the  X  direction 
J  ^ 

F.,  =  body  force  doublet  per  unit  volume  with  same  subscript 

jk 


notation  as  for  T 


jk 


A  =  area 


The  terms  on  the  right-hand  side  represent  the  work  done  by  the  respective 
macro  and  micro  surface  and  body  forces. 

* 

The  internal  energy  U  contained  in  both  macroelements  and  micro¬ 
elements  can  be  expressed  as 


* 

U 


r 

’  J  PMUM 


V 


dv  +  fp'u' 
J 

V 


dV  =  total  energy  -  kinetic  energy 

(5) 


where  ,  U'  are  the  internal  energies  per  unit  mass  for  the 
macroelement  and  microelement ,  respectively. 
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Considering  a  closed  system  consisting  of  both  macroelements  and 
microelements,  thermal  energy  is  transferred  out  of  the  system  at  the 


0  =  J  q  dA 


where  q  is  the  outward  heat  flux.  Conservation  of  mechani  al  energy 


requires 


W  -  Q  =  (T  +  U*) 


The  traction  conditions  at  the  boundary  of  the  macroelement  are: 


tj  =  +  Oy) 


Tjk  +  wijk^ 


where  ni  =  direction  cosines  of  the  outer  normal  to  the  surface 

o ,  .  ,o .  ,  =  components  of  the  respective  symmetric  and  antisym- 
metric  parts  of  the  ordinary  stress  vector  (force/ 
area) 

jj  (i  =  components  of  the  respective  symmetric  and  antisym- 
ijk  ljk 

metric  parts  of  the  doublet  stress  vector  (force 
doublet/area)  acting  normal  to  direction 

Application  of  the  divergence  theorem  and  Eqs.  (8)  yields  the 
following  result  from  Eq.  (4): 


w  -  /  Gy  +  ,i  dV  +  /  [( 

V  L  J  v  L 


(),ijk +  “ijkx'j*  ,i dv 


JVj  dV  +  fr*l 
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where 


[  l.i-M 


i/3X± 


The  heat  transfer  can  be  rewritten  in  terms  of  a  heat  flux  vector 
as  follows: 


Q'  fq)n3  "■  /qj,j  dV  <10) 

A  V 

The  differentiation  on  the  right-hand  side  of  Eq.  (5)  can  be  brought 
under  the  integral  signs  appearing  in  Eqs .  (1)  and  (3)  with  the  following 
results : 


Tt  <T  +  °*>  “  /’Vj  dv  +  I  /*'  dv 

V  V 

+  /(pMVp,i,)  dv 

V 


(11) 


This  interchange  of  the  order  of  differentiation  and  integration  can  be 

6 

justified  by  assuming  small  displacements.  However,,  it  can  also  be 

justified  for  arbitrary  displacements  by  considering  conservation  of 
122 

mass.  Consistent  with  the  assumptions  of  linear  eilastodynamics  the 
densities  are  treated  as  constants. 


The  quantity  u  may  be  written  in  terms  of  its  symmetric  and 

J 

antisymmetric  parts  as  follows: 

"J, i  ■  \  (“l,j  '  ij,i)  '  5  (ii,j  '  "j,!5  =  ‘ij  ■  “lj  (12a> 

where  and  u>^  are  the  macrostrain  and  macrorotation,  respectively, 

given  by  the  following  relations: 


(12b) 
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The  first  two  integrals  appearing  in  Eq.  (15)  vanish  by  virtue  of 
the  following  translational  and  rotational  equations  of  motion: 


°iJ,i+°ij.i+£j-puJ 


(16) 


Uijk,i  + 


pijk,i  +  °jk 


+  F  „  = 


jk 


=  3  p 


dk£^k£ 


(17) 


Thus,  the  only  integral  remaining  in  Eq.  (15)  is  the  last  one,  which 
then  must  be  equal  to  the  right-hand  side,  namely  zero.  Furthermore,  for 
this  integral  to  be  zero  for  arbitrarily  small  regions,  its  integrand 
must  be  zero.  Thus, 


a.  ,t.  . 

ij  ij 


°jkYjk 


+  (p 


ijk 


‘W.i 


-  q 


-  PU 


(18) 


where 


PU  =  MUH+P'U' 


The  microstrain  gradient  is  defined  as  follows: 


Kijk  ^jk,i 


(19) 


Then  Eq.  (18)  can  be  rewritten  in  the  following  form: 


Vu +  V* +  Vijk +  ‘  "j.j  =  ^  <20) 


Equation  (18)  is  the  general  energy  equation  for  an  elastic 
medium  with  microstructure.  When  micros tructurai  effects  are  absent 
(p,  .,  =  u,  .,  =  <k  =  0)  ,  Eq.  (17)  requires  that  the  antisymmetric  part 

of  the  ordinary  stress  tensor  also  vanish  (a.,  5  0)  and  Eq.  (20)  reduces 

5 

to  the  well-known  energy  equation  of  a  classical  elastic  solid. 
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It  is  of  interest  to  note  that  when  all  of  the  stress  (o  and  p)  , 
strain  (e,y)  and  strain  gradient  (k)  components  have  the  same  sign, 
the  presence  of  the  antisymmetric  ordinary  stresses  and  both  the  symmetric 
and  antisymmetric  parts  of  the  doublet  stresses  represent  energy  sinks 
in  the  system. 


It  is  especially  important  to  note  that  even  though  Eq.  (20)  is 
coupled  in  the  sense  that  it  involves  terms  containing  ^ijk*  ^  * 

and  U  ,  there  are  no  cross-energy  terms  such  as  or  ^ijkeij  ’ 

This  considerably  reduces  the  number  of  coefficients  which  are  required 
in  the  constitutive  equations. 


C.  DERIVATION  OF  CONSTITUTIVE  EQUATIONS 

Using  Caratheodory 's  statement  of  the  second  law  of  thermodynamics, 

Boley  and  Weiner^  derived  the  following  relations  for  a  locally  reversible 

* 

thermodynamic  process: 

A  a  +  pTn  =  U  (21a) 

s  s 

Bt6t  -  qj  j  =  pTn  (21b) 

where  ag,  Bt  =  nondissipative  and  dissipative  deformation  variables 

A  ,  B  =  state  functions 
s  t 

p  =  density 

T  =  absolute  temperature 
n  =  entropy  density  function 


Here  we  neglect  viscous  effects  because  it  is  believed  that  shock  wave 
deformation  occurs  so  rapidly  that  there  is  insufficient  time  for  viscous 
forces  to  develop;  thus,  B  =  0  .  Then  Eq.  (21b)  yields: 

q  =  -pTn  (22) 

J  >  J 

*  ”  "  2.23 

According  to  Coleman  and  Mlzel  the  assumption  of  microscopic 

reversibility  is  unnecessary. 


20 


Taking  e^,  j p,  as  state  variables,  we  must  consider  ^ , 


a,.  »  Pj  as  state  functions: 
jk’  ijk 


°ij  "  0ij(£ij’  Yjk’  Kijk,T) 


°jk  ‘=  ajk(Eij’  Yjk’  Kijk>  T) 


pijk  yijk(eij’  V’  Kijk’  T) 


°iieij  +  aijYij  +  PijkKijk  +  pTn  =  pU 


A  free-energy  function  <j>  is  defined  by  the  following  equation: 


**eij’  Yjk*  Kijk’  T)  U(eij*  Yjk*  Kijk’  T)  "  Tn(eij»  Yjk’  T) 


=  U  -  Tn  -  Tn 


But  since 


$  $(eij*  Yjk’  Kijk’  T) 


X  =  li  e  +  — t  y  +  it  k  +  it  f 

*  3e<4  ij  ay.Jjk  3*  *ijk  3T 


Equating  from  Eqs.  (25)  and  (26): 


U  _  Tn  -  Tn  =  it  t  +  — t  v  +  it  ^  +  tt 

U  Tn  rn  ac^ij  +  aYjkYjk  +  3<ijkKijk  +  8T 


Thus,  since  it  has  been  assumed  that  the  other  quantities  are 
independent  of  T  ,  comparison  of  Eqs.  (26)  and  (27)  shows  that 
9  4> 

n  "  ”  9T 
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and 


u - ;t  +  |*  i..  +  |i  +|t  k... 

“ij iJ  3Yjk 3k  SKijk ljk 
Putting  Eq.  (29)  into  Eq.  (23)  gives 


(29) 


+  Vjk  + 


•  /Sdi  • 

^ijk^ijk  =  P^c  £:l 


+  li  j .  + 


ij  iJ  3Yjk  Jk  3Kijk 


"ijk) 


(30) 


Since  it  has  been  assumed  that  <f>  is  independent  of  e^,  Y^j  > 
and  k  ^  »  the  following  additional  equations  must  be  satisfied: 


—  36  *  36  36 

°ij"pae  »°jk-paYjk»  uijk=p  3Kijk 


(31a)  (31b)  (31c) 


Equation  (31a)  has  often  been  used  as  the  thermodynamic  basis  for 
the  symmetry  of  the  elastic  coefficients  in  the  generalized  Hooke's  law 
for  a  general  anisotropic  elastic  solid.  In  analogous  fashion,  Eq.  (31b) 
can  now  be  used  to  establish  the  symmetry  of  the  micros tructural  elastic 
coefficients  in  general  (anisotropic) .  It  is  of  fundamental  importance 
to  note  that  Eqs .  (31)  show  that  thermodynamically  there  can  be  no 
coupling  of  macrostrains  and  microstrains  in  the  constitutive  equations. 
Needless  to  say,  this  simplifies  the  constitutive  equations  considerably. 


To  establish  the  form  of  the  macro  and  micro  constitutive  equations 
more  specifically,  the  free-energy  function  can  be  assumed  to  be  a  poly¬ 
nomial  in  the  macrostrain  components  ,  the  relative  strain  components 

,  the  macrostrain-gradient  components  K-jjk  >  and  the  temperature  T  . 
For  a  linear  material,  defined  as  one  having  linear  constitutive  relations, 
the  highest  degree  polynomial  required  is  of  second  degree.  Thus, 


P<!>  =  a  +  byEy  +  dijkx,eijek£  +  gijeijT 

+b*  Y  .  +  d*  Y..  Y  +  g*  Y..T 

jk  jk  jk  jk  m  &jk  jk 


+B...K..-+D.  „  K,  ,,  K.  +  G  .  K  ,  ..  T 

ij£  ij£  ijk£mn  ijk  £mn  ijk  ijk 


(32) 
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In  view  of  Eqs.  (31)  and  (32),  the  following  constitutive  equations 
can  be  written: 


b .  .  f  d. , , .e , 


+  g^T ;  kl  +  ij 

(33) 

+  g.-.T;  ^  jk 

(34) 

"ijk  "  “ijk  +  DlJkijkKiJk  +  “ijktan^.n  +  G«kT;  *  iJk 


(35) 


There  is  a  total  of  forty  independent  parameters  involved  in  a 
general  three-dimensional  thermoelastic  problem  involving  microstructure. 
These  are  temperature  plus  thirty-nine  mechanical  deformation  parameters: 


Six 
Six 

Twenty-seven  k 


"ij 


(since  »  e^) 
(since  Ykj  =  -Yjk> 


ijk 


(Mindlin8  stated  that  there  are  forty-two  mechanical  deformation  parameters, 
but  he  had  not  determined  that  =  -y^.) 


The  arrays  of  coefficients  have  the  following  sizes: 

(a) 


Coefficient  Array  Type 


b 

d 

g. 


g 

B 

D 

G 


Total  No.  of 
Non-Zero  Coefficients 


No.  of  Independent 
Coefficients 


(a) 

(b) 


IM  6  6 

AM  36  21 

T  3(b)  3(b) 

IM  6  3 

AM  36  15 

T  6  3 

IM  27  27 

AM  729  729 

T  _27  _27 

Total  IM  39  36 

Total  AM  801  765 

Total  T  36  33 

Grand  Total  876  834 

IM  =  initial  mechanical,  AM  =  active  mechanical,  T  =  thermal  expansion 
This  relies  on  the  well-established  fact  that  there  are  no  thermal 
shear  strains  in  linear  theory. 
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The  total  number  of  independent  active  mechanical  coefficients 

listed  in  the  above  tabulation  (765)  is  considerably  smaller  than  the 

10 

1764  claimed  by  Mindlin,  '  although  it  is  still  a  formidable  number. 

The  advantage  of  including  anisotropic  behavior  both  on  a  macroscopic 
and  on  a  microscopic  scale  is  that  it  permits  handling  of  any  material. 
However,  many  metals  and  alloys  of  engineering  importance  (steel,  aluminum, 
some  titanium)  are  macroscopically  isotropic  and  microscopically 
anisotropic  depending  upon  the  crystalline  structure. 

3  24  125 

Recently  Gurtin  ’  proved  a  theorem  stating  that  the  Clausius- 
Duhem  inequality  requires  that  every  elastic  material  be  "simple,"  i.e., 
have  stress,  free  energy,  and  entropy  at  a  material  point  dependent  at 
most  on  deformation  gradient  and  temperature  at  that  point.  However, 
Gurtin  did  not  include  the  possibility  of  any  doublet  stress  effects  in 
his  work.  Thus,  the  present  work  is  not  inconsistent  with  his  theorem; 
rather  it  is  beyond  the  scope  of  his  theorem. 

D.  RELATION  TO  OTHER  THEORIES 

It  may  be  possible  to  relate  the  present  theory  to  the  solid-state 

126-128  128 

theory  of  the  Griineisen  parameter.  For  example,  Gilgarry  has 

127 

argued  that  the  Dugdiile-MacDonald  corrections  for  finite  strain  are 
not  valid;  yet  many  experiments  indicate  better  agreement  with  the 
Dugdale-MacDonald  predictions. 


It  has  been  shown  that  there  is  a  direct  correspondence  between 

129  130 

doublet  stresses  and  a  continuous  distribution  of  dislocations.  * 

This  may  provide  a  basis  for  relating  the  present  work  to  various  theories 

131  132 

of  dislocational  damping. 


The  present  theory  has  no  provisions  for  phase  and  state  transforma¬ 
tions,  which  sometimes  take  place  in  certain  crystalline  materials  under 

3 

strong  shock  conditions. 


24 


III.  DILATION  WAVES  INCLUDING  THE  EFFECTS  OF  MICRO STRUCTURE 


A.  INTRODUCTION 

The  accumulation  of  a  number  of  years  of  experimental  work  with 
shock  waves  in  solids  is  reported  in  Ref.  133.  That  work  has  led  to 
the  conclusion  that  any  proper  explanation  of  shock  wave  phenomena  in 
solids  must  include  the  effects  of  a  rotation  and  deformation  of  microscopic 
regions  even  when  the  input  loading  is  a  longitudinal  pulse.  These 
conclusions  provide  sufficient  reason  for  the  investigation  of  a  theory 
of  "micros tructural"  deformation  in  shock  wave  processes.  A  necessary 
requirement  for  such  an  investigation  is  the  development  of  appropriate 
constitutive  equations  as  done  above. 

These  constitutive  equations  are  used  here  with  the  equations  of 
motion  for  plane  waves. 


B.  DEFINITION  OF  AN  IDEAL  ELASTIC  SOLID  WITH  MICROSTRUCTURE  FOR  SHOCK 
WAVE  STUDIES 

The  definition  of  an  "ideal"  continuum  involves  the  assumption  that 
a  state  variable  (vector)  exists  and  that  a  reversible  path  exists  between 
any  two  possible  conditions  of  the  state  vector.  The  state  vector  for  an 
ideal  elastic  solid  with  microstructure  must  include  the  temperature,  T  , 
and  information  about  either  the  state  of  stress  or  strain  in  both 
conventional  macrostructure  and  the  additional  microstructure.  Constitutive 
equations  have  been  developed  which  involve  the  state  variables  , 

»  anc*  Kijk  (macrostrain,  relative  deformation  between  ma»  ro  and  micro 
media,  and  the  microstrain  gradient),  and  the  state  functions: 


o.ie,,.  y.,  ,  k,  ,  T)  =  symmetric  macrostress  (36) 

ij  ij  jk’  ijk 


o.,  (e  ,  y  i',  >  k4-v»  t)  =  unsymmetric  macrostress  (37) 

lk  ij  ji'  ijk 
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wijk(eij’  Yjk’  Kijk 


>  T)  =  doublet  stress. 


(38) 


This  development  is  valid  whenever  the  kinematics  of  the  physical  problem 
make  Eulerian  and  Lagrangian  coordinates  interchangeable.  In  particulars 
the  assumption  is  made  that 


3u.  3u. 


(39) 


where  u^  =  x^  -  a^  is  the  displacement  of  the  particle,  x_j  is  the 
Eulerian  coordinate  and  a.  is  the  Lagrangian  coordinate.  And  a  similar 
assumption  applies  to  derivatives  in  the  microstructure.  In  the  problems 


where  the  absolute  value  of  displacement  gradients  is  small,  i.e.,  where 


this  is  a  valid  approximation.  However,  we  wish  to  deal  with  shock  waves 
in  the  macrcstructure  where  this  approximation  could  be  considered  invalid. 
The  validity  of  the  assumption  in  reference  to  the  microstructure  also 
needs  to  be  considered.  The  size  of  the  microstructure  elements  has  been 
related  to  the  spacing  (  f  dislocations  in  Ref.  133.  In  Ref.  134,  a 

typical  value  for  the  density  of  dislocations  in  a  crystalline  structure 

9  2  1]  2 

is  suggested  as  10  /cm  before  shock  wave  loading  and  10  '/cm  after 

loading. 


This  would  lead  to  typical  sell  sizes  of  3.10  D  cm/sell  edge  before 
and  3.10  ^  cm/sell  edge  after  shock  wave  loading.  The  same  reference  also 
suggests  that  the  shock  front  thickness  would  typically  be  of  the  order 
of  a  phonon  mean  free  path  or  from  10  ^  cm  to  10  ^  cm,  which  is  comparable 
to  the  sell  size.  This  leaves  the  speculation  that  the  gradient  is  not 
small  in  the  microstructure  with  some  support.  Consequently  it  will  be 
necessary  to  reexamine  the  formulation  of  the  micros tructural 
representation. 


The  Eulerian  equations  of  motion  for  an  elastic  solid  with 
microstructure  are^ 


3^7  <0ij  +  °ij>  +  p£j  '  ’’“j 


(41a) 


3x  (pijk  +  wijk)  +  °jk  +  pFjk  3  p'  ‘Wkfc 


(41b) 


These  equations  are  converted  to  Lagrangian  form  by  substitution  of  the 
Lagrangian  variables. 


9  /-( L)\  :(L)  (L)  _  ;; 

3a±  \  ±i  )  ij  o  j.  Ho  j 

_9 _  /-(L)\  ,  "(L)  :(L)  ,  (L)  _  1  ,  ,  .(L)n2" 

3a  yijkj  Pijk  +  jk  +  P  oFjk  ~  3  P  o(dk)l  )  *1 


(42a) 


(42b) 


The  Eulerian  symmetric  stress  tensor  a,  gives  rise  to  a  non- 


,(L> 


ij 


symmetric  Lagrangian  stress  tensor  .  In  order  to  recapture  a 

symmetric  stress  tensor,  the  Langrangian  tensor  may  be  transformed  to 
the  Kirchhoff  stress  tensor,  and  Eq.  (7a)  becomes 


3  4.  ;oo\  !!il+ ,  f( 

3^'  [v°ik  +  °ikj  3^J+  Pofj 


In  order  to  simplify  the  notation,  we  let  f 


(L) 

=  PoUj 
(L)  _ 


j  -  fj  and  define 


,»)  = 


=  S 


so  that  the  macro  material  equation  of  motion  in  Lagrangian  coordinates 
becomes 


3 

9  a, 


3x 


(Sik  +  Sik)  Saj^ 


+  Pofj  =  Pouj 


(43a) 
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gradient  of  the  scalar 
an  equivalent  "Kirehhoff" 

doublet  stress  in  the  following  form 

(k>  ,  AoN  !i  Ifi 

Pijk  \p  )  3x^  3x^  3x^  Py6e 


ijk 


as  a 


Use  may  be  made  of  the  definition  of  p 
potential  energy  of  an  elastic  solid^  to  obtain 


While  the  relation  between  the  "Kirehhoff"  and  "Lagrangian"  forms  of  the 
doublet  stress  is 


fv .  3x'  3x'  . 

„  00  £  m  _  (L) 

Pi£m  3a^  3a^.  “  Pijk 

The  desired  form  of  the  micro  material  equation  of  motion  is  obtained  as 


3 

3a_. 


7  „  \  3x'  3x^ 

(Mi£m  +  Mi£m) 

'  'Urn 


3x 


+  Sjk  3a^  +  p6Fjk  3  po  (dk£)  \z 


(43b) 


with  the  simplifying  notation  p^  =  M,  d^  =  d  and 


Equations  (43a)  and  (43b)  are  the  equations  of  motion  for  use  with 
plane  wave  phenomena. 


C.  PLANAR  WAVE  IN  MICROMEDIA 

In  a  macromedia  the  planar  wave  is  defined  by  requiring  all  spacial 

derivatives  to  be  zero  except  for  those  in  the  direction  of  the  wave 

3  3 

propagation,  i.e.  t —  (  )  ^  0,  r —  (  )  -  0,  i  =  2,3  .  The  same  idea 

o  d  Si , 

^  3  3 

may  be  extended  to  the  micromedia,  in  which  case  — — -  (  )  £  0,  - — -  (  ) 

dcl^  d3^ 

=  0,  i  =  2,3  .  When  this  is  done,  Eqs .  (43a)  and  (43b)  take  the  form 


(44a) 

(44b) 
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If  we  introduce  the  notion  that  3u/3a  is  a  measure  of  the  wave  strength'^ 
(or  shock  strength),  then  we  may  gi'va  a  meaning  to  3x/3a  in  the  following 

way:  Let  X  =  3u/3a  and  recall  that  u  =  x  -  a  ,  so  that  ~  =  1  -i-  A  , 

then  Eqs.  (44a)  and  (44b)  become  (for  f  =  F  =  0) 


—  [ (S  +  S) (1  +  X) ]  =  pQu 


(45  a) 


o  d 


(1  +  A)S  = 


3  y  3a 


-  4-  [M(l  +  ¥)*] 


(45b) 


Again  following  Bland,  the  stresses  are  replaced  by  derivatives  of  the 

ft 

potential  energy  function,  which  we  will  call  W.  So  that 


=  f-  ^and  e  =  X  +  \  A2^ 


M  = 


3W 

3k 


In  addition,  note  that 


3W  _  de_  _  3W 
3X  3e  3A  3e 


and 


3  h 

lA  32W 

2  2 

3X  31  3  u 

3a  y3 

1 

ro 

3a  3X2  3a2 

then  tq.  (10a)  becomes 


2 

ft  *•  ft  a  •• 

•  ^  +  4-  (U  +  «S]  =  p  u 
3a2  3a 


(46a) 


In  order  to  make  a  similar  elimination  of  M  ,  one  needs  to  differentiate 
the  last  term  in  Eq.  (45b). 

_  -  - 

For  the  present  we  will  not  alter  S. 
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It  is  clear  that  if  Eq.  (46b)  is  introduced  into  Eq.  (46a)  a  very 
complicated  partial  differential  equation  results.  In  order  to  find  any 
acceptable  simplification  of  these  equations,  we  wish  to  know  what 
physical  facts  can  be  used  to  simplify  the  mathematics  of  the  problem  at 
this  point? 


Up  to  the  end  of  the  elaatic  regime,  the  wave  equation  (obtained 
by  s  =  0  in  Eq.  (46a))  gives  a  very  good  approximation  to  the  experimental 
obser  ations.  This  seems  to  be  true  even  when  a  large  number  of 
dislocations  are  present  (i.e.,  when  a  microstructure  is  present). 


30 


We  will  appeal  to  these  observations  as  an  argument  for  setting 
s  =  0  and  see  if  its  consequence  gives  any  meaningful  results.  To  do 
this  we  must  study  the  simultaneous  equations: 


2  2 

3D  3  u  n 

p  u - r  — x  =  0 

°  9X2 


(47a) 


pod* 


V  - 


32W 


a2? 


3k  3a 


_  2  —  —  =  (2  +  ¥)*? 

2  3<  3a  U  ' 


2 
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32W 


3a  3k 


o*  ii  3W 
3a  3k 


(47b) 


and  the  constitutive  equation 

W  =  A  +  BA  +  CX2  +  DAT  +  Fk  +  G<2  +  HkT 


(47c)  * 


with  proper  boundary  and  initial  conditions.  For  now  we  take  D  =  H  =  A  =  0. 
Then  the  behavior  of  Eq.  (47a)  is  well  known  and  Eq .  (47b)  becomes  the 
equation  of  interest.  Substituting  W  into  Eq.  (47b)  leads  to: 


+  2(2  +  404* 


G  +  24- 

3a 


(48a) 


The  right  hand  of  Eq.  (48a)  is  highly  nonlinear.  In  order  to  see  if 
anything  can  be  done  simply,  we  will  arbitrarily  set  all  nonlinear  terms 
to  zero  (the  validity  of  this  step  must  be  evaluated).  So  doing  gives 
Eq.  (48b) 


f  -2G^-|-2f|^=0 
3  -  2  3a 

3a 


(48b) 


Equation  (32). 
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In  order  to  make  the  study  of  Eq.  (48b)  easier,  nondimensionai  variables 
are  introduced  in  the  form: 


2  F  v 


G 

3  =  ? 


d 

C  =7 


with  the  result 


6  -  <J>  -  d  =  0 

44  TT  T 


Equation  (49)  has  the  form  of  the  telegraph  equation  with  space  and 

time  variables  having  their  roles  transposed.  A  substitution  of  the  form 
_ r  /2 

$  =  ^e  transforms  Eq.  (49)  into 

o  =  $  +  y  $  (50) 

44  tt  4 


and  a  formal  solution  with  initial  conditions  of  4>^(4,0)  -  P(4)»  lK4>0) 
=  0  yields 


4>(4»t)  =  e 


J  Jo  (f  V dP 


(51a) 


where 


Q  (p  ,t)  =  |  [P(p  +  t)  +  P(p  -  t)] 


(51b) 


But  P(p  +  t)  corresponds  to  t  <  0  and  does  not  exist.  And  since  the 
initial  velocity  P(4)  is  applied  only  at  the  surface  l  =  0  ,  it  can  be 
represented  as  a  delta  function  P(4)  =  5(4)  so  that  Eqs.  (51)  give 


»(4,t)  =  e~U2  \  J  JQ(j  VT^P2) 


6(p  -  t)  dp 
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The  solution  for  t  >  c,  >  0  is 


♦  «,t>  -{e-r'/2  Jo^i  Vt2  - 


and  for  t  <4  the  signa]  has  not  yet  arrived  at  c,  ,  so  <p  =  0  . 
ing  Eq.  (53)  to  dimensional  notation  yields 


(53) 


Return- 


3u'  1  G  -F/2G  a 

3a'  2  F  6 


0  l2  r  Go  d2 


2  2 

6F  2  F  2 
2  t  2  a 

Gp  d 
o 


(54) 


for  an  initial  microstrain  velocity 


3_  3ul  =  |/  3G3  1 
dt  3a'  f  2p  Fd 


(55) 


corresponding  to  4^  =  1. 


This  result  corresponds  to  an  oscillating  decay  of  microstrain  as 
a  function  of  position  along  the  space  axis. 

If  F  =  0  in  the  constitutive  equation  (Eq.  (47c)),  then  Eq.  (48b) 
becomes 


»od2 

V  -  2G  ~  =  0 
3  a 


.2, 

0 


This  is  a  standard  wave  equation  and  the  microstrain  may  take  any  form 

that  is  a  solution  to  a  wave  equation.  Now  if  either  of  these  cases 

obtains,  then  a  space  oscillation  of  microstrain  is  a  possible  elastic 

solution.  Since  there  has  been  experimental  evidence  of  a  space 

133 

oscillation  in  strain  hardness  after  shock  loading,  it  is  tempting  to 
speculate  on  whether  or  not  these  oscillating  solutions  can  be  related 
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to  the  strain  hardness  observations.  If  the  loading  takes  place  so  rapidly 
chat  the  deformation  extends  its  elastic  range  on  loadings  there  may  be 
some  correlation  between  these  results  and  those  for  which  oscillation  of 
the  strain  hardening  has  been  observed. 


APPENDIX 


THERMODYNAMICALLY  VALID  THERMOELASTIC  CONSTITUTIVE  RELATIONS  CONSIDERING 
TEMPERATURE  DEPENDENCY  OF  THE  ELASTIC  COEFFICIENTS  AND  STRESS  DEPENDENCY 
OF  THE  THERMAL-EXPANSION  COEFFICIENT 

There  are  two  ways  to  formulate  the  thermodynamics  of  solids. 

Traditionally,  it  has  been  assumed  that  deformation  gradients  ^strains) 

and  temperature  are  the  thermodynamic  fluxes  (i.e.,  the  independent 

variables,  also  called  state  variables),  while  stress  was  considered  as 

the  thermodynamic  force  (i.e.,  the  dependent  variable,  also  called  the  state 

19  135 

function).  However,  in  some  instances  *  stress  and  temperature  were 

considered  as  the  independent  variables  and  strain  as  the  dependent 

variable.  In  his  recent  paper  on  thermodynamics  of  strained  solids, 

27 

Kestin  stated  explicitly  that  either  strain  or  stress  could  be  used  as 
the  independent  variable.  The  latter  formulation  is  believed  to  be  mere 
desirable  for  the  following  reasons: 

1.  Prestress  (i.e.,  initial  stress)  is  included  in  a  much  more 

20 

natural  fashion,  as  discussed  by  Burridge  and  Knopoff. 

Prestress  must  be  considered  in  problems  involving  unloading 
waves . 

2.  The  definitions  of  the  elastic  and  thermal-expansion  coeffi¬ 
cients  permit  a  much  simpler  derivation  of  the  interrelationship 
between  temperature  dependency  of  the  elastic  coefficients  and 
stress  dependency  of  the  thermal-expansion  coefficient. 

3.  As  shown  by  Westergaard,  it  is  easier  to  deal  with  nonlinear 
stress-strain  relations  by  complementary  energy  as  a  function 
of  stress  rather  than  strain  energy  as  a  function  of  strain. 

A.  The  yield  criterion  (or  plastic  potential)  in  plasticity  theory 
is  a  function  of  stress  not  strain. 
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In  view  of  these  considerations,  especially  the  second  one,  the 
formulation  having  stress  and  temperature  as  the  independent  variables 
will  be  used  here. 


The  energy  equation  is  as  follows: 

ei/ij  +  T*  =  ”  (Al) 

where  e..  is  the  strain  tensor,  a.,  is  the  stress  tensor,  T  is 

ij 

absolute  temperature,  n  is  the  entropy  per  unit-reference-state  volume, 

* 

U  is  the  internal  energy  per  unit-reference-state  volume,  and  a  dot 
denotes  a  derivative  with  respect  to  time. 

The  free-energy  function  is  now  defined  as  follows: 


Ho±y  T)  =  U(o..,  T)  -  Tn(a..,  T) 
Putting  Eq.  (A2)  into  Eq.  (Al)  gives 


(A2) 


(e..  -  3*/3o  )o  -  (n  +  3*/3T)T  =  0  (A3) 

6  124  125 

Since  thermodynamically  *  <}> ,  n  ,  and  z .  .  are  independent 

.  x3 


of  0 .  .  and  T , 
ij 


=  3-ti/9o^  ,  n  =  — 3  4>/9T 


(A4)  (A5) 


Onsager's  reciprocal  relation  implies  sufficient  continuity  of  <f> 
that  the  order  of  partial  differentiation  can  be  interchanged.  Thus,  the 
following  equality  must  be  satisfied: 

■k 

Since  the  Lagrangian  formulation  is  used,  this  unit- reference-state 
volume  basis  is  equivalent  to  using  a  unit-mass  basis,  but  it  has  the 
advantage  of  eliminating  the  symbol  for  refere->ce-state  density  in 
numerous  subsequent  equations. 
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0 


ij 


or 


3 

ifii 

3 

!!m  | 

3T 

T  =  3T 

,8u  1 

(A6) 


To  obtain  constitutive  equations  which  are  linear  in  stress,  the 
free-energy  function  is  assumed  to  have  the  following  form: 


9  =  B1(T)J1  +  B2(T)J2  +  C1(J1,T)(J1)(T  -  1q)  +  F(T) 


(A7) 


where  Tq  is  the  temperature  at  the  relaxation  state  (i.e.,  the  tempera¬ 
ture  at  which  the  material  is  free  of  strain) ,  and  B2  are 
temperature-dependent  elastic  coefficients,  is  the  usual  thermal- 

expansion  coefficient,  F  is  an  arbitrary  function  of  temperature  only 
(detailed  form  selected  to  represent  the  nonlinear  dependence  of  specific 
heat  at  constant  volume  on  temperature)  and  and  J2  are  stress 

invariants  defined  as  follows: 


where  i, 
given  by: 


J1  -  «ii*  J2  * 


J® 


(A8) 


j,  i,  m  =  1,  2,  3  and  6^  is  the  generalized  Kronecker  delta 


+1  if  i,j  form  an  even  permutation  of  £,m 
-1  if  i,j  form  an  odd  permutation  of  £,m 
0  if  i,j  do  not  form  a  permutation  of  £,m 


The  form  of  the  C..  term  in  Eq .  (A7)  implies  the  existence  of  a 

1  25 

relaxation  state,  a  point  which  has  been  questioned  by  Eckart,  who 

offered  no  specific  alternative  to  use  in  constructing  constitutive 

relations . 
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In  view  of  Eq.  (A4) 


£jLj  =  0*/3Jx)  (3J1/3o  )  +  (3*/3J2)(3J2/3o  )  (A9) 


Differentiation  of  Eq.  (A8)  yields 


3J./3o.  .  =  6 . . ,  6J_/3o .  .  =  6  .J,  -  o. . 

1  ij  ij  2  ij  ij  1  ij 


(A10) 


where  6  is  the  simple  Kronecker  delta  (takes  values  of  unity  for 
i  =  j  ,  zero  for  i  #  j). 


Inserting  Eqs.  (A7)  and  (A10)  into  Eq.  (A9) ,  one  obtained 


ciJ  -  “l(T)Jl4ij  +  B2<n<Vl  ‘  V 


+  t. .(T  -  T  )3(C.J.)/3J1 
IJ  o  1  1  1 


(All) 


or 


en  =  2B1(T)o11  +  [2BX(T)  +  B2(T)](o22  +  o^) 


+  (T  -  Tq)3(C1J1)/3o11 


(A12) 


(Similar  relations  for  e22  and  e^)  • 


and 


e12  “  "B2(T)al2 


(A13) 


(similar  relations  for  z^  and  e^) 


Comparing  Eqs.  (A12)  and  (A13)  with  those  of  classical  elasticity 
theory,  we  find 
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2B^  -*■  E  ,  where  E  is  Young's  modulus 
2B^  +  B^  -v/E  ,  where  v  is  Poisson's  ratio 
3(CiJi)/3J^  -*■  a  ,  lineal  coefficient  of  thermal  expansion 

To  check  whether  the  Onsager  requirement,  represented  by  Eq.  (A6) 
is  satisifed,  the  mixed  partial  derivatives  are  taken  as  follows: 


d  B  d  B 

2  IT  6ij  h  Tt  (<$ij 


=  i- 

3T\30i  J 


1}  +  6ij  t(T  "  Va(ClJl)/3Jl] 


(A14) 


Equation  (A14)  shows  that  Eq.  (A6)  is  satisfied  without  placing  any 
restrictions  whatsoever  on  the  form  of  the  constitutive  functions 
B1(T),  B2(T),  C^.T)  ,  and  F(T) . 

The  specific  entropy  is  computed  from  Eqs.  (A5)  and  (A7)  with  the 
following  result: 


n  ■»  -J:  d  B./d  T  -  J.  d  B_/d  T  -  J.  3  [  (T  -  T  )C.]/3T  -  d  F/d  T 
11  2  2  1  o  1 


(A15) 


Since  there  can  be  no  entropy  production  during  continuous  deformation 
of  a  perfectly  elastic  material,  the  specific  entropy  must  be  independent 
of  the  stresses.  Thus,  from  Eq.  (A15) : 


d  B2/d  T  =  0 


(A16) 


J1  d  B2/d  T  +  3 [ (T  -  To)C1]/3T  =  0 


(A17) 


It  can  be  shown  that  B2  is  the  reciprocal  of  the  shear  modulus 
G  =  E/2(l  +  v)  ;  thus,  Eq.  (A16)  implies  that  the  shear  modulus  should  be 
independent  of  temperature.  However,  it  is  an  experimentally  observed 
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fact  that  the  shear  modulus  _is.  temperature  dependent.  This  can  be  accom¬ 
modated  quite  easily  by  adding  a  term  of  the  following  form  to  the  free- 
energy  function  <j>  in  Eq.  (A7). 

C2(J2,T)(J2)(T  -  Tq) 

In  the  stress-strain  equation,  Eq.  (All),  this  adds  a  term  of  the  form 
<V  -  ‘’ijX'1  -  T0)3(C2J2)/3J2 

which  is  analogous  to  the  deviatoric-thermal-expansion  effect  considered 
theoretically  and  experimentally  by  Dillon. ^  Addition  of  this  term  does 
not  prevent  satisfaction  of  Eq.  (A6). 

To  obtain  some  useful  information  from  Eq.  (A17) ,  use  is  made  of  the 
following  relation  observ  d  in  experiments  by  Rosenfield  and  Averbach: 

Vv  T)  -  C10(I)  -  cii'T>0kk 

Rosenfield  and  Averbach  used  only  uniaxal  loading  and  measured  C|2  only 
in  the  direction  of  that  loadirg .  However,  it  seems  intuitively  apparent 
that  a  stress  in  direction  1  could  affect  the  thermal  expansion  in  direc¬ 
tion  1  only,  i.e.,  since  the  loading  is  directional,  its  effect  on 
should  also  be  directional.  This  probably  could  be  proved  by  functional 
analysis.  However,  the  present  theory  is  an  isotropic  one;  thus,  the 
preceding  equation  must  be  rewritten  as  its  isotropic  equivalent: 

Cl(Jl't)  ■  C10(C)  +  Cll(t)Jl  <A18) 

Then  Equation  (A17)  leads  to 

Jx  dBL/dT  Td [ (T  -  To)C1Q]/dT 

+J1  d[(T  -  To)cn]/dT  =  0  (A19) 
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Again  invoking  independence  of  stress,  Eq.  (A19)  implies: 


d[(T  -  To)C10]/dT  =  0 


(A20) 


and 


dBn/dT  +  d[(T  -  T  )C11]/dT  =  0  (A21) 

i  o  ii 

Integration  of  Eq.  (A20)  yields: 

(T  -  T  )C10(T)  =  Const.  (A22) 

o  iu 

Equation  (A22)  does  not  appear  to  be  realistic;  this  suggests  that  further 
work  on  this  topic  is  required. 

Equation  (A21)  expresses  an  interrelationship  between  and  . 

Since  the  experiments  of  Rosenfield  and  Averbach  were  conducted  at  room 
temperature  only,  the  validity  of  Eq.  (A21)  cannot  be  checked  with 
experimental  data  at  present.  Table  1  lists  some  sets  of  specific  forms 
for  and  B^  which  satisfy  Eq.  (A21). 

Since  B^  =  (2E)  ,  Eq.  (2)  in  the  body  of  the  report  can  be  put 

in  the  following  form: 

Cn  =  2  dB1/dT  (A23) 

Thus,  Eq.  (A23)  corresponds  to  a  g  value  of  2,  where  g  is 
defined  as  follows: 

g  =  Cn/(dB1/dt)  (A2A) 

It  is  noted  that  the  form  of  Eq.  (A23)  is  identical  to  Cases  2-A  in 
Table  1,  even  though  the  coefficient  g  differs  considerably.  It  is 
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TABLE  1 


SOME  SETS  OF  SPECIFIC  FORMS  OF  COEFFICIENTS  B 

EQUATION  (A21) 


AND  C  WHICH  SATISFY 


Set  No. 
1 

2 

3 

4 


"11 

K,/(T  -  T  ) 
l  o 


K.  -  K0T 
4  3 


K 


K5e 


3 

-k/T 


k6t 


g (T)  =  C11/(dB1/dT) 


-1 

-[1  +  (T  -  T  )k/T2]-1 
o 

-[1  +  (T  -  To)p/T]-1 


•k 

K^,  k,  p  are  constants 

noted  that  for  temperatures  near  to  the  strain-free  temperature  (tQ)  , 
Cases  2-4  predict  g  =  -1  .  However,  Rosenfield  and  Averbach's  experiments 
for  five  metals  gave  a  value  of  g  within  3%  of  2  (although  individual 
variations  for  specific  metals  differed  from  2  by  6%  to  100%). 

Apparently  Rosenfield  and  Averbach  made  a  mathematical  error  in 
applying  Eq.  (A6) ,  and  in  doing  so,  they  obtained  the  interrelationship 
Eq.  (2).  This  is  contrary  to  current  thinking  in  thermodynamics  that 
Onsager's  relation  should  not  place  numerical  restrictions  on  the 
constitutive  equation  coefficients,  but  rather  should  be  used  only  in 
determining  the  general  form  of  the  constitutive  equations.  Nevertheless, 
it  is  puzzling  that  Eq.  (A23)  agrees  with  experimental  results  better  than 
with  the  values  »,f  Cases  2-4  in  Table  1. 
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